A Hot DQ White Dwarf in the Open Star Cluster M35 
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ABSTRACT 

We report the discovery of a hot DQ white dwarf, NGC 2168:LAWDS 28, that is a 
hkely member of the 150- Myr old cluster NGC 2168 (Messier 35). The spectrum of the 
white dwarf is dominated by C II features. The effective temperature is difficult to esti- 
mate but likely > 20 000 K based on the temperatures of hot DQs with similar spectra. 
NGC2168:LAWDS 28 provides further evidence that hot DQs may be the "missing" 
high-mass helium-atmosphere white dwarfs. Based on published studies, we find that 
the DBA WD LP 475-242 is likely a member of the Hyades open cluster, as often as- 
sumed. These two white dwarfs are the entire sample of known He-atmosphere white 
dwarfs in open clusters with turnoff masses > 2Mq. Based on the number of known 
cluster DA white dwarfs and a redetermination of the H-atmosphere:He-atmosphere 
ratio, commonly known as the DA:DB ratio, we re-examine the hypothesis that the H- 
to He-atmosphere ratio in open clusters is the same as the ratio in the field. Under 
this hypothesis, we calculate that five He-atmosphere WDs are expected to have been 
discovered, with a probability of finding fewer than three He-atmosphere white dwarfs 
of 0.08, or at the f« 2a level. 

Subject headings: white dwarfs — open clusters and associations: individual (NGC 2168) 



1. Introduction 

White dwarfs (WDs) have one of two types 
of atmospheres: those dominated by hydro- 
gen (spectral class DA), and those dominated 
by helium (the non-DA spectral classes). The 
ratio of the numbers of DA to non-DA stars 
is a function of temperature, culminating in 
the "DB gap", ranging from ^ 45 000 K to 
PS 30 000 K (e.g., Liebert et al. 1986). Un- 
til recently, no He-atmosphere WDs were 
known to exist in this temperature range, 
though several have now been discovered in 
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the Sloan Digital Sky Survey (Eisenstein et 
al. 2006b). The temperature dependence of 
the DA to non-DA ratio suggests that some 
He-atmosphere WDs must change their atmo- 
spheric composition, an effect often ascribed 
to the diffusion of hydrogen to the surfaces of 
cooling DO stars and the later mixing of this 
hydrogen veneer into the thicker He-envelope 
due to the appearance of an atmospheric con- 
vection zone (Fontaine & Wesemael 1987). 

Kalirai et al. (2005) note that even though 
a few dozen WDs are known to exist in open 
star clusters with turnoff masses greater than 
2Mq, no DB WDs are known in these clus- 
ters. However, this study excludes other 
WD spectral types with He-dominated atmo- 
spheres, such as the DBA WD LP 475-242 



(WD 0437+138), a potential member of the 
Hyades. Kalirai et al. (2005) also give only 
a brief mention of the impact of the DB gap 
without calculating the impact of the chang- 
ing DA to non-DA ratio in this temperature 
range. Subsequent to their publication, anal- 
ysis of the SDSS DR4 WD sample (Eisenstein 
et al. 2006b) has found that the DA to non- 
DA ratio is larger than the value of 3.5 used in 
Kalirai et al. (2005), warranting a re-analysis 
of the cluster non-DA deficiency. 

In this Letter, we report the discovery of 
a hot DQ (He-dominated atmosphere with 
opacity dominated by atomic carbon) WD in 
the field of the young open cluster NGC 2168 
(Messier 35). Several DA WDs in this clus- 
ter have already been identified and character- 
ized (Wilhams, Bolte, & Koester 2004). The 
identification of this He-atmosphere WD pro- 
vides insight into both the origin of DQ WDs 
and the apparent deficiency of He-atmosphere 
WDs in open clusters. 

2. Observations 

We obtained UBV imaging of the open 
cluster NGC 2168 on UT 2004 January 27 
with the Mosaic- 1 camera on the Kitt Peak 
Mayall 4-m telescope. Photometry was de- 
rived via point-spread function fitting using 
the DAOPHOT II program (Stetson 1987). 
We select objects with blue excess as candi- 
date WDs as discussed in Williams, Bolte, 
&: Koester (2004). The color-magnitude di- 
agram for confirmed cluster WDs is shown in 
Figure 1. NGC 2168 has an age of 150 ± 60 
Myr (von Hippel 2005) and [Fe/H] = -0.21 ± 
0.10 (Barrado y Navascues, Deliyannis, & 
Stauffer 2001). The main-sequence turnoff 
mass is ~ 4.5M0 (Girardi et al. 2002); based 
on the initial-final mass relation, WDs in 
NGC 2168 are expected to be massive (e.g., 
Weidemann 2000; Ferrario et al. 2005). 

The open diamond in Figure 1 is NGC 



2168:LAWDS 28 (hereafter LAWDS 28), 
a(J2000) = 6'^08'"13';50, 5(J2000) = +24°20'32'.'5, 
and we assign it the McCook & Sion (1999)- 
style designation WD J0608+2420. This 
object has a F-magnitude of 21.63 ± 0.03, 
B-V = 0.02±0.04, and U-B = -1.10±0.04. 
As its photometric properties (open diamond 
in Figure 1) are similar to those of known 
WDs in M35, we selected this object for 
follow-up spectroscopy. 

The spectroscopic data were obtained on 
UT 2005 November 26 using the blue camera 
of the Low-Resolution Imaging Spectrometry 
on the Keck 1 telescope (Oke et al. 1995). We 
placed a l" slit on the object at the paral- 
lactic angle and selected the 400 line mm^^, 
3400A blaze grism to obtain a spectral reso- 
lution of ~ 6A. Our total exposure time was 
4800s, split into four 1200s exposures. The 
images were combined and the spectrum was 
extracted and calibrated (wavelength and rel- 
ative flux) using standard IRAF routines. 

The extracted spectrum is peculiar and 
markedly different from other spectra ob- 
tained in our open cluster WD survey, both 
for confirmed WDs (DA, DB and DCs) and 
the non-WD sample contaminants (mostly 
QSOs, hot subdwarfs, and field horizontal 
branch stars). The spectral features are con- 
sistent with C II spectral features, leading us 
to conclude that this object is a DQ WD. 

The spectrum for LAWDS 28 is shown in 
Figure 2 (top), along with the SDSS 2.5- 
m spectrum for SDSS J010647.93+151327.9 
(bottom), the object most similar from the 
sample of hot DQ stars showing C I, C II, O I 
and/or O II multiplets in Liebert et al. (2003). 
The tick marks label the approximate posi- 
tions of multiplets of C II taken from Moore 
(1972). Included are the weaker multiplets 
omitted from the plots in Liebert et al. (2003). 
Both multiplets with strong and weak labo- 
ratory intensities appear in LAWDS 28; the 
4267A line, often observed to be strong in hot 



stars, is present in LAWDS 28 but is not the 
only dominant feature. While a case could be 
made that the SDSS object shows additional 
features, possibly due to O II in its spectrum, 
the cluster WD is matched almost perfectly 
with C II transitions alone. It is possible that 
the strong C III multiplet near 4647-51 A may 
contribute weakly, but no C I or oxygen lines 
are matched except where they would blend 
with the C II. 

Hot DQ WDs are thought to have atmo- 
spheres dominated by helium. Thejll et al. 
(1990) analyze the atmospheric abundances of 
the hot DQ G35-26 (WD 0203+207) and find 
a carbon abundance nc ~ 0.03 and a hydro- 
gen abundance of riu < 0.01, implying that 
nue ~ 0.96 (Liebert et al. 2003). The hot 
DQ G227-5 (WD 1727+560) also was found 
to have only trace abundances of carbon and 
hydrogen (Wegner & Koester 1985). There- 
fore it is likely that LAWDS 28 has a He- 
dominated atmosphere despite showing no ev- 
idence for helium features in the spectrum. 
The lack of such features in many hot DQs 
is thought to be due to the masking of he- 
lium features by the dominant carbon opacity 
(Liebert et al. 2003). 

3. Discussion 

3.1. Cluster membership, mass, and 
temperature 

Given the lack of a model fit to the spec- 
trum of LAWDS 28, it is difficult to determine 
many fundamental parameters of this WD 
with precision. However, the available obser- 
vations provide a preponderance of evidence 
that LAWDS 28 is a hot, massive cluster- 
member WD. We now outline this evidence. 

The distance and reddening of NGC 2168 
are somewhat uncertain. Kalirai et al. (2003) 
adopt E{B — V) = 0.20 from Sarrazine et al. 
(2000) to derive (m - M)o = 9.80 ± 0.16 
from BV photometry, while Sung &: Bessell 



(1999) calculate (m - M)o = 9.60 ±0.1 and 
E{B-V) = 0.255 ± 0.024 from UBVRI pho- 
tometry. Either case leads to an apparent dis- 
tance modulus (m — M)y w 10.4. At this 
distance, the absolute magnitude of LAWDS 
28 is My = 11.2; assuming the Sarrazine et 
al. (2000) reddening value, {B-V)o = -0.22 
and (U — B)q = —1.23. This absolute mag- 
nitude is consistent with evolutionary models 
of DA and DB WDs with Mwd > O.8M0 and 
hotter than 20 000K (Bergeron, Wesemael, & 
Beauchamp 1995). However, DQ photometric 
indices are likely to differ from those of DAs 
and DBs, and in fact the hot DQs in Liebert 
et al. (2003) are bluer in u — g than pure- 
He atmosphere models. As seen in Figure 1, 
LAWDS 28 is also bluer than the given DB 
models. 

The effective temperature, T^s, of LAWDS 
28 cannot be determined accurately, as no 
spectral models for these hot DQs exist. 
The B — V color of this WD is bluer than 
any DB model from Bergeron, Wesemael, &: 
Beauchamp (1995), though these models are 
not calculated for T^s > 30 000 K. By match- 
ing the observed B — V color to reddened 
blackbody models calculated using the syn- 
phot package in STSDAS (Bushouse & Simon 
1994), we derive a blackbody temperature of 
22 200t^^[5[] K. For the hot DQs in Liebert 
et al. (2003), SDSS colors were converted to 
a "pure helium" effective temperature, Tne, 
based on colors for pure helium-atmosphere 
WDs with log 5 = 8. For SDSS J0106+1513, 
the comparison star shown in Figure 2, Thc = 
20 000 K, while the DQs SDSS J0005-1002 
and SDSS J0236-0734 have Tue = 29 000 and 
30 000 K, respectively (Liebert et al. 2003). 
Both of these stars show C I in addition to 
C II in their spectra; the lack of C I in the 
spectrum of LAWDS 28 would argue that 
this DQ is hotter than the other two. While 
the continuum slope of SDSS J0106+1513 is 
bluer than that of LAWDS 28, this does not 



necessarily mean that LAWDS 28 is cooler, 
as the spectra have not been corrected for 
extinction. Since neither the blackbody tem- 
perature nor the pure-helium temperatures 
take the spectral features of the DQs into 
account, these are only educated guesses as 
to Tefj. Given the very blue colors and the 
ion ratios, we surmise that LAWDS 28 is a 
hot (Teff > 20 000 K) WD. We also note that 
the coolest DA WDs in NGC 2168 (similar in 
B-V to LAWDS 28) have T^s ~ 30 000 K. 

The discovery of a hot DQ white dwarf in 
this young cluster is consistent with the case 
reiterated in Liebert et al. (2003) that such 
WDs may be massive. Thejll et al. (1990) 
found a high gravity for the hot DQ star G 
35-26 from their detailed atmospheric analysis 
and fit to the spectral lines of carbon, hydro- 
gen and helium. This has been confirmed by 
a trigonometric parallax measurement (Dahn 
et al. 2006). A parallax is also available from 
the same source for the similar field star G 
227-5, the spectrum of which was analyzed 
by Wegner & Koester (1985), and it also ap- 
pears to have a mass near 1 Mq. It is ex- 
pected from theoretical evolutionary calcula- 
tions of asymptotic giant branch stars that the 
helium-shell burning consumes a greater frac- 
tion of the helium layer with increasing mass 
(Kawai, Saio, & Nomoto 1988), with further 
calculations leading to a similar conclusion in 
Thejll et al. (1990). A smaller helium enve- 
lope as the star enters the white dwarf cooling 
phase in turn means that dredge-up of car- 
bon (and perhaps oxygen) from the edge of 
the core into the convective envelope and at- 
mosphere should occur at a higher tempera- 
ture - in the DB temperature range, in which 
LAWDS 28 apparently is. If all hot DQ stars 
are in fact massive, this may account for the 
absence of massive DB stars in the sample of 
Beauchamp (1995, see also Beauchamp et al. 
1996); this is the only large sample of field DB 
stars for which detailed determinations of the 



relevant parameters (Tefj, log 5, and therefore 
mass) are available. Nearly all of these DB 
stars have masses near the canonical peak of 

0.6 Mq. 

Finally, hot DQs are very rare. In the SDSS 
DR4 WD catalog (Eisenstein et al. 2006a), 
only 0.3% (28 out of 9316 white dwarfs) are 
hot DQs; with a surface density of approxi- 
mately one for every 180 square degrees. We 
estimate the number of hot white dwarfs in 
our 0.5 sq. degree imaging of NGC 2168 by 
integrating the WD luminosity function from 
the Palomar-Green Survey (Liebert, Berg- 
eron, & Holberg 2005) for My > 11 to a 
volume corresponding to limiting apparent 
magnitude of V = 21.6. This estimate there- 
fore allows for the possibility that LAWDS 
28 is less massive than IMq, and therefore 
brighter and background to NGC 2168. The 
estimated number of intrinsically bright WDs 
in our imaging area is found to be ~ 10. If 
0.3% of these are hot DQs, we would expect 
0.03 hot DQs in our imaging area. Thus, the 
likelihood that LAWDS 28 is an interloping 
field DQ in our 0.5 square degree field of view 
is extremely low. However, until a distance 
modulus or proper motion can be determined 
for LAWDS 28, some uncertainty about its 
membership will remain. 

In summary, the photometric and spec- 
troscopic data for LAWDS 28 are consistent 
with it being a hot {T^fi > 20 000 K), massive 
(M > O.8M0) WD located at the distance of 
NGC 2168. Until an accurate model atmo- 
sphere can be fit to the data or until proper 
motion information is available, it is not pos- 
sible to state conclusively that LAWDS 28 is 
a member of NGC 2168. 

3.2. Other non-DA WDs in open clus- 
ters 

LP 475-242. — LP 475-242 (WD 0437+138) 
is a DBA white dwarf about one core ra- 
dius (3 pc) from the Hyades cluster cen- 



ter. Beauchamp (1995) fits its spectrum 
to get Tcff = 15 355 K and log 5 = 8.26, 
suggesting a modestly above-average mass 
(0.751 ibO.OSSMo) for a DB white dwarf. Hy- 
drogen lines were included to get a trace hy- 
drogen abundance of log(nH/?^He) = —4.5, so 
the atmosphere is clearly He-dominated. This 
WD has been used in some cluster data com- 
pilations to determine the initial-final mass 
relation, most recently that of Ferrario et al. 
(2005). 

LP 475-242 was hsted by Luyten (1971) as 
a possible white dwarf member of the Hyades 
on the basis of its proper motion, faint mag- 
nitude, and blue color. Greenstein (1976) 
determined a photometric parallax of 0'.'023, 
in close agreement with the Hyades distance. 
Eggen (1985) and Weidemann et al. (1992) 
consider LP 475-242 to be an established 
Hyades member. Accurate proper motion 
measurements (Table 1) agree almost exactly, 
in size and direction, with motion expected for 
a Hyades member. Finally, under the assump- 
tion that LP 475-242 is at the systemic veloc- 
ity of the Hyades, the resulting gravitational 
redshift leads to a mass of 0.748 ± O.O37M0 
(Reid 1996), consistent with the spectroscopic 
mass (Beauchamp 1995) and with the masses 
of DAs in the Hyades (Bergeron, Liebert, & 
Fulbright 1995). It therefore seems highly 
likely that LP 475-242 is a Hyades member. 

LB 3600. — It is also possible that a 
third cluster DB is already known, LB 3600 
in Messier 67 (Fleming et al. 1997). Based 
on their spectral fits, this WD is about one 
magnitude too bright to be a member of 
M67, though their spectrum is very noisy, 
and DB spectra are difficult to fit with pre- 
cision. As the deficiency of open cluster DBs 
noted by Kalirai et al. (2005) was specifically 
for clusters with main-sequence turnoff masses 
> 2Mq and the M67 turnoff mass is Ki IAMq, 
we do not consider this object in our analysis 
below. 



3.3. The open cluster DA:non-DA ra- 
tio 

Kalirai et al. (2005) note the lack of DB 
WDs in open clusters. Though there are no 
known "pure" DB WDs in open clusters, this 
overlooks the physical reality that non-DA 
spectral types in the DB temperature range 
- DBs, DBAs, DBZs, and hot DQs - have at- 
mospheres dominated by He. As He has a very 
low opacity, trace contaminants in the atmo- 
sphere can easily dominate the opacity, mask- 
ing He features. In order to study whether 
the ratio of hydrogen-dominated atmosphere 
WDs to helium-dominated atmosphere WDs 
differs from star clusters to the field, these 
other spectral types should be included in any 
analysis. Given our discovery of LAWDS 28, 
a likely He-dominated atmosphere WD in an 
open cluster, given the probable Hyades DBA 
LP 475-242, and given a new determination of 
the field DA:DB ratio from Sloan Digital Sky 
Survey data (Eisenstein et al. 2006b) , we now 
revisit the issue of the H- to He-atmosphere 
ratio in open clusters. 

A new measurement of the field DA:DB 
ratio. — The large numbers of field WDs 
discovered during the Sloan Digital Sky Sur- 
vey has permitted a new determination of the 
DA:DB ratio (including DBAs) (Eisenstein et 
al. 2006b). For Teg < 20 000 K, DA:DB« 
5, while for 20 000 K < Teg < 45 000 K, 
DA:DB~ 12.5. Because of the near degener- 
acy of DB spectra over the temperature range 
of 20 - 28 000 K, it is difficult to determine 
the lower T^s boundary of the deficit, though 
it is most likely in the 20 — 22 000 K range. 
Thus, the "DB-gap" is really a deficiency, not 
a complete lack, of WDs in the temperature 
range 45 000 K < Tcff < 20 000 K, somewhat 
larger than previous estimates of the DB gap 
range. As this is a large, uniform sample of 
WDs, and as the DA:DB ratio is provided as 
a function of temperature, we adopt these ra- 
tios for calculating the number of expected 



He-atmosphere WDs in the open cluster sam- 
ple. 

Expected numbers of He-atmosphere WDs 
in the open cluster sample. — For the sam- 
ple of open cluster WDs, we use the sample 
as presented by Ferrario et al. (2005). This 
work applies stringent membership criteria, 
so few non-cluster members should contam- 
inate the sample. We calculate the number 
of non-DA WDs using the Eisenstein et al. 
(2006b) DA:DB ratio in two temperature bins: 
Teff < 20 000 K (DA:DB= 5; 21 open cluster 
DAs) and 20 000 < T^s < 45 000 K (DA:DB= 
12.5; 17 DAs). We thus expect ^ 5.6 non-DA 
WDs in the current open cluster WD sam- 
ple. This number assumes that the Eisenstein 
et al. (2006b) DA:DB ratio would be unaf- 
fected by inclusion of non-DB and non-DBA 
He-atmosphere WDs in the calculations. 

Given two probable He-atmosphere cluster 
WDs (LAWDS 28 and LP 475-242), we use 
Poisson statistics to calculate the likelihood 
of finding two or fewer non-DA WDs to be 
0.08. If only one of these He-atmosphere WDs 
is a cluster member, this probability drops to 
0.02. This is akin to 1.9a and 2.3a, respec- 
tively. Thus, while the ratio of DA to non-DA 
WDs in open clusters is low, it is not greatly 
inconsistent with the expected number. With 
coming larger samples of open cluster WDs, 
the question of whether the DA to non-DA 
ratio in open clusters is the same as that in 
the field can be addressed with much more 
certainty. 
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Table 1 
Proper motion measurements of LP 475-242. 



fJ-a 



l^s 



^ 



PA Reference 



Comments 



0.093^ -0.036^ 0.10 111 (1) 

0.100 -0.021 0.102^ 102^ (2) 

0.097 -0.023 0.100^ 103'^ (3) 

0.094 -0.018 0.096^ 101^ (4) 

0.098^" -0.021^ 0.100^ 102^ (5) 



LP 475-242 = Reid 400 

LP 475-242 = NOMAD 1039-0046673 

Predicted motion for a Hyades member at position of LP 475-242 



'^Calculated from data in reference 



Note. — Proper motions are in arcseconds per year; position angles are degrees east of north. 



References. — (1) Luyten 1971, (2) Luyten, Hill, & Morris 1981, (3) Reid 1992, (4) Zacharias et al. 2004, 
(5) Perryman et al. 1998 




Fig. 1. — Color-magnitude diagram of con- 
firmed white dwarfs in NGC 2168. Curves 
are evolutionary tracks for pure H- (i.e., DA; 
solid) and pure He-atmosphere (i.e., DB; 
dashed) white dwarfs hotter than 10 000 K 
with masses of 0.6Mq, 0.8Mq, I.OMq, and 
1.2M0(top to bottom, respectively). Tracks 
have been moved to the cluster distance and 
reddening. Filled squares are the confirmed 
cluster member white dwarfs from Williams, 
Bolte, &: Koester (2004). The open diamond 
with la error bars is NGC 2168:LAWDS 28; 
the error in V is smaller than the point size. 
This white dwarf has photometry consistent 
with known cluster white dwarfs and with 
massive white dwarf evolutionary models. 
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Fig. 2.— Spectrum of NGC 2168:LAWDS 28 
(top) compared to the known hot DQ SDSS 
J010647.9-M51327.9 (bottom). Dashed tick 
marks indicate the location of C II features. 
The qualitative match of the two spectra and 
the fact that most or all of the features in 
the top spectrum match the C II features lead 
us to conclude that NGC2168:LAWDS 28 is a 
hot DQ white dwarf. 



